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r ^ 1. Abstract 

^ H We study stability of occurrence of black holes and naked singularities that arise as a final 

state for a complete gravitational collapse of type I matter field in a spherically symmetric N 
dimensional spacetime with equation of state p = kp, < k < 1. We prove that for a regular 
initial data comprising of pressure (or density ) profiles at an initial surface t — ti, from 
which the collapse evolves, there exists a large class of the velocity functions and classes of 
^ . solutions of Einstein equations, such that the spacetime evolution goes to a final state which 

is either a black hole or a naked singularity. We further prove that in an infinite dimensional 
separable Banach space, the set of regular initial data leading the collapse to a black hole or 
\Q ' a naked singularity, forms an open subset of the set of all regular initial data. In this sense, 

the gravitational collapse leading to either a black hole or a naked singularity is stable. These 
results are discussed and analyzed in the light of the cosmic censorship hypothesis in black 
CN ■ hole physics. 
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1 Introduction 



The arrival of Tolman's dust solutions pQ in 1934 to Einstein's field equations in spherically sym- 
metric spacetime has led to investigation of nature of singularity in various models like dust, null 
dust, perfect fluid, quark matter etc. in later years. Considering Tolman solution, Oppenheimer 
and Snyder [2], in 1939, have studied the collapse of spherically symmetric homogeneous dust 
cloud, which led to the general concept of trapped surfaces and the black hole (BH). In 1951, 
Vaidya identified the general metric appropriate to a spherically symmetric distribution of null 
dust [3] and considering linear mass function, Papapetrou and Kurada [4] have shown that the 
central singularity will be a naked singularity and persistent, when the collapse is sufficiently 
slow. Hereafter, many numerical and theoretical studies were carried out in different models by 
various authors [5] presenting counter examples to Cosmic Censorship Cenjecture(CCC). 

In a recent paper, R. Goswami and P. S. Joshi [6] pointed out that the occurrence of naked 
singularity(NS) as end state of dust collapse beginning from smooth initial data with physically 
relevant restrictions, can be removed when one goes to higher dimensional space-time. This has 
restored CCC in dust collapse for sufficiently higher dimension of the space-time. Many other 
cases of the dust collapse have been studied by various authors [7] . 
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Lot of importance is being given to the physically relevant models such as Type I matter fields 
that include most of the known physical forms of matter like dust, perfect fluid, etc. The authors 
in [8] have discussed the role of initial data in spherically symmetric gravitational collapse for 
Type I matter fields. 

Joshi and Goswami have studied the case of an isentropic perfect fluid with linear equation 
of state in four dimensional spacetime, wherein, they pointed out that the occurrence of NS/BH 
evolving from regular initial data depends on the choice of rest of the free functions available 
[9j, and also they have proved that in N dimensional dust collapse, the occurrence of NS can be 
removed when one goes to higher dimensional space-time, thus restoring CCC |10| . 

In recent years, special attention is given to higher dimensional gravitational collapse since 
the development in string theory indicate that gravity may be a higher dimensional interaction 
[llj . Hence, it would be interesting to know the status of CCC in higher dimensional spacetime. 
These studies may help in possible appropriate mathematical formulation of the conjecture or 
assist in providing any possible proof of CCC [12]. N. Dadhich et.al. have studied spherically 
symmetric collapse of fluid with non- vanishing radial pressure in higher dimensional space-time 
with equation of state p r = kp through construction of the root equation governing the nature 
(BH versus NS) of the central singularity [13J. 

However, all of the initial data is not independent, so it poses a question, which initial data 
set shall evolve into a NS as a result of gravitational collapse? Further, having a NS for certain 
initial data set is not enough because if it looses its characteristic for small perturbation in the 
neighbourhood of that data set then the NS is not serious enough to challenge cosmic censorship 
conjecture |14| . Therefore, the all important question is that, "is the NS developed from certain 
initial data set a stable one?" . We study this behaviour of NS using definitions of stability and 
genericity following Abraham and Marsden [15\ [16] . 

Introduction of equation of state has importance as it makes the collapsing models physically 
more relevant. Does the linear equation of state p = kp stimulates the outcome of NS/BH phases 
as end state of gravitational collapse? We investigate this question through analysis of field 
equations of N-dimensional gravitational collapse with linear equation of state in section 2. The 
study of radial null geodesies in TV dimensions is carried out in section 3. In section 4, we analyze 
the condition of existence of naked singularity to prove our assertion mentioned above by using 
existence theory of first order ordinary differential equations [17] . Thereafter, in section 5, we 
prove openness of initial data set and then we discuss genericity of this set in section 5.1 



2 Collapse with linear equation of state 



The general spherically symmetric metric describing N dimensional space-time geometry of a 
collapsing cloud can be described in comoving coordinates 
(t,r,9\9 2 ,..,e N - 2 ) by 



ds 2 = - e Mt,r) dt 2 + e mt,r) dr 2 + R 2^ r)dn 
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is the metric on an [N — 2) sphere. The stress-energy tensor for Type I field in a diagonal form 
is given by [18] 



-p, T r r = T 0l = p . 



(2) 



The quantities p and p are the energy density and pressure respectively. We take the matter field 
to satisfy weak energy condition which implies p>0; p + p>0. Linear equation of state for perfect 
fluid is 

p(t, r) = k p(t, r) where k 6 [0, 1] . (3) 



Einstein field equations for the metric ([T]) are derived as 

_ (N - 2)F' _ 1 (N - 2)F 
P ~ 2R N ~ 2 R' ~~k 2R N ~ 2 R 



(4) 



,'=- — [ln(p)f (5) 
R'G - 2RG u' = (6) 

G - H = 1 -Jr-z- (7) 

where F = F(t, r) is an arbitrary function, and has an interpretation of the mass function for 
the cloud. It gives the total mass in a shell of comoving radius r on any spacelike hypersurface 
t = const. The energy conditions impose condition on F that F(t, r) > 0. In order to preserve 
the regularity at the initial epoch, we have F(ti,0) = 0, that is, the mass function should vanish 
at the center of the cloud. Since we are considering collapse, we have R < 0, i.e. the physical 
radius R of the collapsing cloud keeps decreasing in time and ultimately reaches R = 0, which 
denotes the singularity where all matter shells collapse to a zero physical radius. The functions 
G and H are defined as G(t,r) = e~ 2 ^R' 2 and H(t,r) = e~ 2u R 2 . 

Let us define a new function v(t, r) by v(t, r) = R/r, and now we use the scaling independence 
of the comoving coordinate r to write 

R(t, r) = r v(t, r) 

where v(ti,r) = 1; v(t s (r),r) = 0; v < (8) 

The time t = t s (r) corresponds to the shell- focusing singularity at R = 0, where all the matter 
shells collapse to a vanishing physical radius. The introduction of the parameter v as above allows 
us to distinguish the spacetime singularity from the regular center, with v = 1 at the initial epoch, 
including the center r = 0, which then decreases monotonically with time as collapse progresses 
to value v = at the singularity R = 0. At the regular center the mass function F(t, r) behaves 
suitably so that the density remains finite and regular there at all times till the occurrence of 
singular epoch. A general mass function for the cloud can be considered as 

F(t,r) =r N ~ 1 M(r,v) (9) 

where M(r, v) is regular and continuously twice differentiable. Using equation ([9]) in equation 
dH), we obtain 

N- 2 (N — 1)M + r[M >r + M tV v'} (N - 2)M„ 

P ~ 2 X v N - 2 (v + rv') ~ 2k v N - 2 ' [ } 

Then as v — > , p — > oo and p — > oo i.e. both the density and pressure blow up at the singularity. 
At the initial epoch, the regular density distribution takes the form 

N -2 

Po(r) = —z—[rM(r,l),r + {N-l)M(r,l)] ■ (11) 



We rearrange equation (jlUj) as follows 

krM ir + Z(r,v)M tV = -(N-l)kM (12) 

where 

Z(r,v) = (k + l)rv' + v . (13) 

Equation (|12p has many classes of solutions but only those classes of solutions should be consid- 
ered which satisfy energy conditions, which are regular and gives p — > oo as v — > 0. This means 
the energy conditions and equation of state p = kp isolate the class of functions M(r, v) so that 
the mass function F(t,r), the metric function u(t,r) and the function b(r) (to follow) evolve 
as the collapse begins according to the Type I field equations. The energy condition imposes 
constraint on M(r, v) asM„ < 0. Hence we set 

M ^v) = J^ (14) 

where M Q is a positive constant, such that energy density p — > oo as v — > 0. On integrating 
equation ([5]), we obtain the general metric function, 

u(r,v) = -l[ln(p)]. (15) 

where I = ^xj. Define a suitably differentiable function A(r,v), A{r,v) t v = u'/R' so that 
equation (fl~5j) takes the form 

A (r,v), v = -^j- (16) 



lMr,v), v ] v=1 = --^ (17) 



At the initial epoch, we have 

j P» 
Po{r) 

whereas using equation (|10p . the relation between the function M and A, at all epochs is given 
by 

(JV - 2)M,, 



A(r,v) v Rl = -I 



In 



2k v 



N-2 



From above, it is clear that M(r, v) is a decreasing function with respect to v and if we consider 
a smooth initial profile for the density in such way that density gradient vanishes at the center 
then we have A(r, v) = rg Q (r,v) where g (j",v) is another suitably differentiable function. 
Also, the use of definition of A(r, v) in equation ([6]) yields 

G{t,r) = d(r) e 2rA (19) 

where d(r) is another arbitrary continuously differentiable function of r. Following comparison 
with the dust model, we can write 

d(r) = 1 + r 2 b(r) (20) 

where b(r) is the energy distribution function for the collapsing shells. Define a function h(r, v) 

as 

h (r,v) = ^^- (21) 
and substituting this equation, together with (fT9|) and ([20]) in equation fl7J), we get 



(N-3)/2 i = _ p -l v N-3 h ^ + bu JV-3 e 2rA + M ( 2 2) 
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where negative sign is chosen since, for the collapse, v < 0. Integrating the above equation , we 
have 

t(v,r)= / = dV . (23) 

In above equation, the variable r is treated as a constant. Expanding t(v,r) around the center, 
we get 

t(v,r) = t(v,0) +r X (v) + 0{r 2 ) (24) 



where the function 

() _dt_ 1 r 1 y( N - a V 2 B t r(p,v) 

X{V} -d^ r=0 - 2 J v 0(0, v) 3 / 2 



dv (25) 



and 

B(r, v) = p~ 21 [v N - 3 h(r, v) + bv N ~ 3 e 2rA + M(r, v)] . (26) 
The time when the central singularity develops is given by 

...f^. (27) 



/o B(0,«)i/2 

The time for other collapsing shells to reach the singularity can be expressed by 

t 8 (r) = t(0,r) = t So + r X (0) + 0(r 2 ) . (28) 
where x(0) = nm x( v )- (29) 

Now, it is clearly seen that the value of x(0) depends on the free functions b(0), M(Q, v) and 
h(0, v), which in turn, depend on the initial data at the initial surface t = U. Thus, a tangent 
to the singularity curve t = t s (r) is completely determined by the given set of density, pressure, 
velocity function v and function b. Further, from equation (|22p . onu = const, surface, we can 
write 

v (N - 3)/2 v' = x{v)B{0,v) 1/2 + 0(r) . (30) 

Let us analyze the nature of x(0) from the point of view of equation (|25p . The existence 
of the limit depends on the cumulative effect of all the terms present in the equation but since 
0(0, v) > as v — > 0, provided M(0,0) > and v € [0, 1], therefore positive or negative sign of 
x(0) absolutely depends on the sign of B >r (0, v). 

Now, B(0, v)>0 and if 

)r .(O, v) = _ 2 W*) g (M +p -2i [ v N-z {h ^ v) +b / (0)) + M, r (0,t;)] < 

then integrand in equation (|25l ) is positive and hence, x(0) > 0. This is possible, since density 
is very high in the neighbourhood of the singularity and it decreases as r increases from r = 0, 
therefore, p'(0,t) should be negative to be physically reasonable and hence sign of 0, r (O, v) 
depends on the sign of [?j Af_3 [h ) r (0, v) + b'(0)] +M r (0,u)] i.e. depends solely on the initial 
data. The introduction of an equation of state has ensured that the system of Einstein equations 
is closed and so freedom to specify any function is exhausted. Therefore, the functions h(r, v) be 
written in terms of b(r) and M(r,v). Hence, the choice of class of functions b(r) and M(r,v) 
shall decide the final fate of gravitational collapse. So, the occurrence of NS or BH phases exist, 
depending on the choice of the class of function M(r,v). 



3 Radial null geodesies 



Now, we investigate below when there will be families of null geodesies emanating, which will be 
future directed and outgoing, and which terminate in the past at the singularity, thus making 
the communication from the singularity to an outside observer possible, as opposed to a black 
hole situation. Hence, for examination of the nature of central singularity at R = 0, r = 0, we 
consider the equation of outgoing radial null geodesies, given by, 

dt 
dr 

Further, we write the null geodesic equation in terms of the variables ( u = r^,R), choosing 

l(N-l) 



(3 = rrr[ l^^ — k], and using equation 0, we have 




(N-3)-(N-l)k 

v v 2 

[(iV-3)-(iV-l)fc] 

a i 2 




(31) 



If the null geodesies terminate at the singularity in the past with a definite tangent, then at the 
singularity, we have dR/du > 0, in the (u, R) plane with a finite value. All other points r > on 
the singularity curve are covered since F/R N ~ 3 — > oo with dR/du — > — oo and only the singularity 
at the centre r = could be naked. For the case, when R' > near the central singularity, we 
have 

, , R dR 
x n = lim lim 



t-»t s r->o u du 
and use of equation (f30|) for v', equation (j3Tj) yield 



(32) 

t — >t s ,r— Yo 



k 
h+l 



(iV-Dd-fc) (l-k)(N-l) r— ( 2 . 

and the radial null geodesic emerging from the singularity in (R, u) co-ordinates is R = x Q u, or 
in (t, r) plane, it is given by 

t-t s {0)=x o r 13 . 

Therefore, x Q > iff x(0) > 0, and hence x(0) > is a sufficient condition for the occurrence 
of the NS at the center of the cloud as the end state of gravitational collapse of a sufficiently 
dense star when it loses its equilibrium state, and continual collapse begins with regular initial 
data of density and pressure profiles p = kp. There will be radially future outgoing null geodesies 
emanating from the singularity, giving rise to a locally NS at the center. However, if x(0) < 
then we have a black hole solution, as there will be no such trajectories coming out. Also, it is 
observed from above equation that parameter k ( because (1 — k) > 0) and dimensionality TV" 
do not take part in the decision making of the sign of x(0). Thus the final fate of collapse of a 
massive star leads to either a BH or a NS with appropriate choice of initial data for N > 4. For 
x(0) = , we will have to take into account the next higher order non-zero term in the singularity 
curve equation, and a similar analysis can be carried out by choosing a different value of f3. 

4 Existence of function b(r) leading to NS 

In this section, we prove the assertion mentioned in the introduction. We choose b(r) to satisfy 
the differential equation on a constant u-surface 

l v (N-3)/2 B M) 

2 g(r ,„)3/2 =*fo») 



a 



for < r < r c where B(r, v) is a continuous function defined on D such that B(0, v ) < for all 
v in [0, 1]. It will then follow that 



x(o) 



lim x( v ) 



B(0,v)dv > 0. 



(34) 



This condition ensures that central shell-focusing singularity will be naked. 

We, now, discuss the existence of b(r) as a solution of the differential equation (1331) , We put 



x(r, v) = p~ 2 1 [v N ~ 3 h(r, v) + bu iV ~ 3 (l + r z h(r, v)) + M(r, v) 
a continuous function of r, in equation (|33p and rearranging it, we get 



.N-3, 



dx 
dr 



1 



v (N-3)/2 



2B{r,v)x 3/2 \ = f{x,r) 



with the initial condition 

x(0, v) = p(0, v)~ 2 1 [v N ~ 3 h(0, v) + b(0)^- 3 + M(0, v)} . 



(35) 



(36) 



(37) 



Let us ensure the existence of C 1 - function x(r, v) as a solution of above initial value problem 
defined throughout the cloud. The function f(x, r) is continuous in r, with x restricted to a 
bounded domain. With such domain of r and x, f(x, r) is also C 1 - function in x which means 
f(x, r) is Lipschitz continuous in x. Therefore, the differential equation (I36p has a unique solution 
satisfying initial condition (j37|) . 

Further, we can ensure that the solution will be defined over the entire cloud i.e. for all r in 
[0,r c ] by using the freedom in the choice of arbitrary function B(r,v). For this , we consider the 
domain [0, r c ] x [0, d] for some finite d. Let S be the supremum of the modulus of f(x,r). Then 
the differential equation (j36|) has a unique solution defined over the entire cloud provided [T7] 



<mf(r c ,-) = -. 



This yields 



max 

0<r<r c ,0<x<d 



1 



2B(r, v)x 



3/2 



v (N-3)/2 

Condition (|38p . in turn, will be satisfied if the weaker condition 



d 

< —. 



< \B(r,v)\ x 3/2 < 



2r c 



(38) 



(39) 



(40) 



holds for all r in [0, r c ]. 

The collapsing cloud may start with r c small enough so that the expression 
which is always positive, satisfies the condition (|40p with x restricted to a bounded domain. 
We then have infinitely many choices for the function B(r,v), which is continuous and satisfies 
conditions (|40|) and B(0, v) < for each choice of v. For each such B(r, v), there will be a unique 
solution x(r, v) of the differential equation (j36j) . satisfying initial condition ([37|) . defined over the 
entire cloud and in turn, there exist a unique function b(r) for each such choice of B(r,v), that 
is given by the expression 



b(r) 



v N - 3 (l + r 2 h(r,v)) 



x{r, v) 



N 



2 / 



(N - l)M + r(M r + M )V v') 
v N ~ 2 (v + rv') 



2 l 



,N-3 



h(r, v ) — M(r, v) 
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over [0, r c ]. Thus, for a given constant v -surface and given initial data of mass function F(t,r) = 
r N ~^M{r, v) and h(r, v) = (e r s i2 ( r ' 1 ') — l)/r 2 satisfying physically reasonable conditions (expressed 
on M), there exists infinitely many choices for the function b(r) such that condition B(0, v) < 
is satisfied. The condition continues to hold as v — > 0, because of continuity. Hence, the central 
singularity developed in the collapse is a naked singularity. 

5 Stability of Naked Singularity 

The analysis in section 2 shows that the only conditions on the initial data which evolves the 
collapse into a naked singularity are the energy conditions. Hence, for stability, we have to 
examine energy conditions only. What we now show is that, the set of such initial data satisfying 
energy conditions forms an open subset of the space of all initial data. For this , we use the 
technique of Saraykar and Ghate [20J. The energy conditions are p>0 and p + p>0 but since 
p = kp, p>0 is the only sufficient condition required in this set up and can be written in the form 

[(N - 1)M + rM r + rv'M )V ] =E 1 >0, 

on the domain D. Here, we have used an assumption that v is an increasing function of r so that 
v ' > on any surface t = tj say, so that no shell-crossing singularity condition (ie.B! > 0) holds. 
So, this is a valid assumption. 

We assume that X be an infinite dimensional Banach space of all C 1 real-valued functions 
defined on T>, endowed with the norms 

|| M(r,v) ||= sup|M| +sup|M jr | +sup|M„| 

V V ' V 

These norms are equivalent to the standard C 1 norm 

|| M(r,v) ||= sup(|M| + |Af r | + |M >W |) (41) 
v 

Let G = {M(r, v) : M > 0, M is C 1 , E x > on V} be a subset of X. 

We show that Q is an open subset of X. For M in Q, let us put 5 = min(M), 7 = min(i^i), and 
A3 = max(u') for varying r in [0,r c ] and v e[0, 1], the functions involved herein are all continuous 
functions defined on a compact domain T> and hence, their maxima and minima exist. We define 
a positive real number 

1 ■ fx 7 7 7 i 
U = — mini 0, —, r, , — k 

^ 2 1 ' 3(N - 1)' 3r c ' 3r c \ 3 r 

Let Mi(r,v) be C l in V with \ M — M\ ||< /i. Using definition pi]) , we get \M\ - M\ < p, 
\Mi >r — M r \ < p and \M\ tV — M >v \ < p over D. Therefore, for choice of p, the respective 
inequalities are 

Mi > M - - > 0, (42) 

(N — l)\Mi — M\ < — , (43) 
6 

r\M lr - M r \ <r c \Mi r - M r \ <-, (44) 

6 

rv'\M 1;V - M ;V \ < r c X 3 \M hv - M c | < 1 (45) 



that are satisfied on T>. Combining inequations (|43I44I45|) from above, we have 



(N - 1)| Mi - M| + r|M 1>r - M r | + rv'\M 1>v - M tV \ < ~ < 7 < E 1 

Therefore, | [(N - 1)M 1 + rM 1)T . + rw'Mi^] - ^ | < #i where #1 > on V. 
Hence, [(N - l)Mi + rMi, r + ru'M 1)V ] > on V. 

Thus, Mi > 0, Mi is C 1 and [(iV - l)Mi + rM hr + ru'Mi,,,] > throughout V. Therefore, 
Mi(r, v) also lies in Q and hence, Q is an open subset of X. Now, let Y be an infinite dimensional 
Banach space of all C 1 real-valued functions defined on T>. The set H of functions h(r, i>) would 
form (by similar technique) an open subset of Y. Then Q x % is an open subset in the product 
space X xY. 

Taking (Mi, hi(r, v)) in the neighbourhood of (M, h(r, u)) inGx?i, and using equation (|40p 
analogously for Mi and hi, we have a choice of infinitely many B\(r, v), such that for each such 
-E>i(r, v), there will exist a unique bi(r) so that the initial data of mass function r 3 Mi and hi(r, v) 
together will lead the collapse to formation of a naked singularity. Thus, naked singularity arising 
from (M, h(r, v)) is C 1 - stable in the sense defined above but since in earlier analysis the initial 
data functions are at least C 2 , therefore NS arising from initial data is C 2 -stable. 

The analysis given above, guarantee the existence of a metric function u(t, r) for a given 
initial data set. Such choice of v{t, r) and expressions for G and H together will yield the metric 
([1]) as an exact solution leading to the occurrence of naked singularity. Thus, given initial data 
and the weak energy conditions above determine the matter distribution and the metric of the 
space-time. 

5.1 Issue of Genericity of Black Holes and Naked Singularities 

We have seen that in Type I matter field collapse, the end state of collapse is governed by the 
choice of initial data. We have also shown that the set of initial data leading the collapse to a 
naked singularity or a black hole forms an open subset of the full set of initial data, and thus 
both the occurrences are stable. 

It is also clear that these sets as subsets of initial data set leading the collapse to BH and NS 
are mutually disjoint. Since both of them are open, none of them can be dense in X x Y. Thus, 
in the strict Mathematical sense |15l 116] . both these occurrences are not generic. 

Nevertheless, each of these sets are substantially big sets. Following the measure theoretic 
analysis on infinite dimensional spaces given in Section IV of [21], it can be argued that these 
sets have a non-zero measure. 

6 Conclusions 

Our main conclusions are as follows: 

1(a). Given a C 2 - mass function M(r,v) and a C 1 - function h(r,v), on any v = constant surface, 
we find a C 1 - energy distribution function b(r) such that the collapse ends in a naked singularity. 
1(b). With physically reasonable conditions put on mass function M, and with C 1 - function 
h(r, v), the initial data consisting of (M, h(r, v)) leading the collapse to a naked singularity forms 
an open subset of the space X x Y. This establishes stability of naked singularity with respect 
to initial data. 

2. Similar analysis for the case x(0) < ( that is B(0,v) > for all v in [0, 1]) shows that the 
occurrence of black hole is stable but not generic. Thus, both the occurrences , namely, black 
holes and naked singularities, are stable but not generic. 



n 



3. Also, it is found that equation of state parameter k and dimensionality N of spacetime do not 
take part in the decision making of the sign of x(0). So for critical value x(0) > 0, a spectrum of 
radial null geodesies can emanate from the central singularity but many of them may fall back 
to the singularity. Further, it can be investigated whether k and N assist in making a locally NS 
to be a globally NS which can be seen by a faraway observer. 
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